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We employ 100 ns molecular dynamics simulations to study the influence of cholesterol on structural and dynamic 
properties of dipalmitoylphosphatidylcholine (DPPC) bilayers in the fluid phase. The effects of the cholesterol 
content on the bilayer structure are considered by varying the cholesterol concentration between and 50 %. We 
concentrate on the free area in the membrane and investigate quantities that are likely to be affected by changes in 
the free area and free volume properties. It is found that cholesterol has a strong impact on the free area properties 
of the bilayer The changes in the amount of free area are shown to be intimately related to alterations in molecular 
packing, ordering of phospholipid tails, and compressibility. Further, the behavior of the lateral diffusion of both 
DPPC and cholesterol molecules with an increasing amount of cholesterol can in part be understood in terms of 
free area. Summarizing, our results highlight the central role of free area in comprehending the structural and 
dynamic properties of membranes containing cholesterol. 



I. INTRODUCTION 

Cholesterol is one of the most prominent molecular species 
in the plasma membranes of mammalian cells. It is a tremen- 
dously important molecule, a component essential f or the 
very existence and multiplication of cells jFinegoldl Il993t 
Iphvo-Rekila et al., 2002, and references therein). It is abun- 
dant in the plasma membranes of higher organisms: depend- 
ing on the exact lipid composition, the plasma membr ane may 
conta in of the order of 20 - 50 % cholesterol ( Alb erts et all 
11994 . 

Eukaryotic cells do not seem to be able to grow and dif- 
ferentiate properly without cholesterol. It has been firmly 
established that cholesterol mod ulates the physical proper- 
ties of the plasma membrane ( McMuUen and McEIhanevl 
[1996). A finite cholesterol content has been said to im- 
prove the characteristics of a simple phospholipid bilayer 
and allow for wider variations in th e lipid composition of 
the membrane llVist and Davislll990h . Perhaps not surpris- 
ingly, cholesterol is one of the primary molecules in lipid 
rafts ( Edidin, 2003: Silvius, 200.3; Simons and lkonen. 1997> 
and references therein), i.e., microdomains rich in cholesterol, 
sphingomyelin, and saturated phospholipids. Rafts have been 
thought to confine proteins involved in e.g. signal transduction 
events, and hence act as platforms for adhesion and signaling. 
Consequently, one could well imagine that as cholesterol al- 



ters the properties of the bi layer, it might affect the func tioning 
of the embedded proteins dCantoj. 1 1 999t I Yeaelel [1991 . 

The effects of cholesterol on the properties of phospho- 
lipid bilayers are diverse. In the physiologically relevant fluid 
phase, adding cholesterol to the bilayer leads t o increased 
orientational order in the phospholipid tails ("Chiu et all 
|20p2; HofsaBetal., 2003; McMullen and McElhanev, 1991 
ISankara m and Thomcson, 1990b) and smaller average areas 
per molecule (.Petrache et aU .1999) . In other words, choles- 
terol modifies the packing of molecules in bilayers. Other im- 
portant effects are changes in passive permeability of small 
solutes (.Jedlovszkv and Mezei. .2003; Xiang. .1999. and ref- 
erences therein) and suppressed later al diffusion of phos- 
pholipids in bilayers with cholesterol ( Almeida et al., 1992; 
[ Gallaetal. , 1979; HofsaBetal., 2003; Poison etal., 2001; 
IVattulainen and Mouritsentl2003l) . Both permeability and lat- 
eral diffusion, in turn, are strongly affected by the amount and 
distribution of free volume or area in a membrane, i. e., space 
not occupied by phospholipids, cholesterols, or water. Choles- 
terol thus seems to simultaneously influence packing, free 
area, diffusion, and permeability in lipid bilayers, and it is 
reasonable to expect that the changes in these properties are 
somehow coupled. 

Although there is a wealth of information on the effects 
of cholesterol on lipid bilayers, the interplay of packing, free 
area, diffusion, and permeability has not yet been studied sys- 
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tematically. Experimental electron density profiles ("Mcintosh', 
17978) and deu terium nuclear magnet ic resonance (NMR) 
data (Sankara m and Thompson! Il990bl) suggest that choles- 
terol should influence the packing inside membranes. Fluo- 
rescence recovery after photobleaching (FRAP) experiments, 
in turn, have been used to study the dependence of late ral dif- 
fusion coefficients on free area ( I Almeida et all fT99^ . More 
information at the atomic level, however, is essential for gain- 
ing a detailed understanding of the effect of cholesterol on 
lipid bilayers. Such atomic-level information can be obtained 
from computer simulations. Molecular dynamics in particu- 
lar provides a unique tool to investigate both the structure and 
dynamics of lipid membranes with a level of detail missing 
in any experimental technique. Until recently, however, sys- 
tematic simulation studies have been limited by the extensive 
computational requirements. 

In the present study, we investigate the cholesterol- 
induced changes in packing, free area, ordering, and lat- 
eral diffusion in phospholipid bilayers. Specifically, we 
study the presumptive interplay between these changes. To 
this end, we employ 100 ns molecular dynamics simula- 
tions on dipalmitoyl phosphatidylcholine (DPPC) / cholesterol 
bilayers, with cholesterol concentrations ranging from 
to 50 mol %. While detailed multi-nanosecond simula- 
tion stu dies on atom ic level have emerged only very re- 
cently jHofsaB et all 120031: IScotll l2002t iTieleman et^ 
Il997h . there exist large amounts of experimental studies 
for DPPC / cholesterol bilayers j McMullen and McElhanev, 



122^; ISankaram and Thompson! Il990albt: iVist and Davisl 
1990t and references therein). These previous studies and the 



experimental results in particular offer us an excellent plat- 
form for comparison. 

In order to further enhance the understanding of the effect 
of cholesterol on bilayers, we introduce a novel method for 
investigating the packing and free area in bilayers. The scope 
of this technique is very wide. It allows us to estimate how 
much space DPPC, cholesterol, and water molecules on aver- 
age occupy in different regions of the bilayer. Consequently, 
it yields information on the amount and location of free space 
in the bilayer As discussed below, this is related to various 
structural aspects such as the ordering of lipids in a membane. 
Our method also provides valuable insight into dynamic prop- 
erties. For example, our approach allows us to determine the 
area compressibility across a membrane, and hence yields in- 
formation on rate-limiting regions for lateral diffusion. In ad- 
dition, as the method enables us to examine changes in free 
area with an increasing cholesterol content, we may estimate 
diffusion coefficients in terms of free area theories for lat- 
eral diffusion. The present approach can be applied to a wide 
range of different kinds of membrane systems, including one- 
and multi-component bilayers, and bilayers with embedded 
solutes, probes, and proteins. 

We find that cholesterol strongly affects the amount of 
space occupied by molecules in different parts of a phos- 
pholipid bilayer The close-packed areas occupied by the 
tails of DPPC molecules can be explained by the orde ring 
of the tails, and a simple relation ( Pet rache et all Il999h can 
be used for quantifying the dependence of close-packed area 



on ordering. The amount and location of free space is sig- 
nificantly reduced by an increasing cholesterol content, and 
clearly reflect the total space occupied by DPPC and choles- 
terol molecules. The lateral diffusion coefficients, too, show a 
substantial decrease with an increasing choles terol concentra- 
tion. We find that so-called free area theories (lAkneida et all 
119921 ICohen and Turnbulll Il959^ iGalla et all Il979h .' which 
are essentially two-dimensional mean-field models, correctly 
predict this reduction, but are not applicable to quantitatively 
describing lateral diffusion in lipid bilayers. 



II. MODEL AND SIMULATION DETAILS 

We studied fully hydrated lipid bilayer systems consisting of 128 
molecules, i. e., DPPC's and cholesterols, and 3655 water molecules. 
Since the main focus of this paper is on studying the effects of choles- 
terol on phospholipid bilayers, we were interested in bilayers with 
varying amounts of cholesterol. To this end, we studied a pure DPPC 
bilayer and composite DPPC / cholesterol bilayers with six different 
cholesterol molai' fractions: x = %, 4.7 %, 12.5 %, 20.3 %, 29.7 %, 
and 50.0 %. 

The starting point was a united atom model for a fu lly hydrated 
pure DPPC bilayer that has been validated previously IPatra et all 
2003 : Tieleman and Berendsen. 1996). The parameters for bonded 
and non-bonded interactions for DPPC molecul es wer e taken from 
a study of a pure DPPC bilayer tBerger et^l Il997h available at 
|http : / /mo ose .bio.ucalgary. ca/Downloads/lip id. itp| 
The partial charges are from the underlying model descrip- 
tion (.Tieleman a nd Berendsen. 1996) and can be found at 
[http : / /moose ■ bio . ucalqar y . ca/Downl oads/dppc . i tp| 
For water, the SPC model iBerendsen et aU .1981) was 
used. As our initial configuration for the pure DPPC 
bilayer we used the final structure of run E discussed 
in Ref. (Tieleman and Berendsen, 199^ and available at 
http://moose.bio.ucalgary. ca/Downloads/dppcl 2 8 ■ pdb| 
The bilayer is aligned such that it lies in the xy-plane, i.e., the 
bilayer normal is parallel to the z-axis. 

The cholesterol force field and the initial shape of 
an individual cholesterol molecule were taken from 
http : / /www . gromacs .org/topologies/ uploadedjnolecu 
cholesterol . tgz (Holtie et al., 2001). Cholesterols were in- 
troduced to the bilayer by choosing DPPC molecules from the pure 
phospholipid bilayer at random and replacing them by cholesterols. 
The same number of DPPC molecules was replaced in each of the 
two monolayers. In practice, the center of mass (CM) of a choles- 
terol molecule was moved to the CM position of the removed DPPC 
molecule. The main axis of inertia of each inserted cholesterol was 
parallel to the z-axis, and each molecule was rotated by a random 
angle around the z-axis. 

The molecular dynamics (MD) simulations were perform ed at a 
temperature T = 323 K using the GROMACS I Lindahl et all 1200 ll) 
molecular simulation package. The time step for the simulations was 
chosen to be 2.0 fs. The lengths of all bonds were kept constant with 
the LINCS algorithm iHess et allll997l) . Lennard- Jones interactions 
were cut off at 1.0 nm without shift or switch functions. Long- 
range e lectrostatic interactio ns were handled using the Particle-Mesh 
Ewald ( 'Essman et all 1 199^ method, which has been shown to be a 
reliable metho d to account for l ong-range interactions in lipid bi- 
layer systems iPatra et all l2003h . Electrostatic interactions within 
1.0 nm were calculated at each time step, while interactions beyond 
this range were determined every t en time steps. These choices fol- 
low the parametrization of DPPC iTieleman and Berendseri Il99d) 
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and correspond to a scheme called twin-range cutoff. 

After an inital energy minimization, we needed to equilibrate the 
system to fill the small voids left by replacing DPPC molecules 
by somewhat smaller cholesterol molecules. The equilibration was 
commenced by 50 ps of NVT molecular dynamics with a Langevin 
thermostat using a coupling time of 0.1 ps, i.e., every 0.1 ps the ve- 
locities of all particles were completely randomized from a Maxwell 
distribution corresponding to the target temperature. This complete 
loss of memory after 0.1 ps reduces the amount of ballistic motion 
of atoms inside a void. The equilibration was continued by 500 ps 
of NpT molecular dynamics at a pressure of 1 bar with a Langevin 
thermostat and a Berendsen barostat (Berendsen et al., 198^. The 
time constant for the latter was set to 1 ps, and the height of the sim- 
ulation box was allowed to vary separately from the cross-sectional 
area of the box. 

Finally, for every cholesterol concentration, we performed 100 ns 
of MD in the NpT ensemble with a Berendsen thermostat and baro- 
stat iBerendsen et all 1X9841 The barostat was the same as the one 
described above, and the thermostat was set to separately couple the 
DPPC, cholesterol, and water molecules to a heat bath with a cou- 
pling time of 0.1 ps. The six simulations took a total of approx- 
imately 60,000 hours of CPU time. For all systems up to and in- 
cluding the cholesterol molar fraction of 29.7 %, a simulation time of 
100 ns guarantees a good sampling of the phase space. The results for 
50 % cholesterol should be regarded with some caution, as the dif- 
fusion of the DPPC and cholesterol molecules is already quite slow, 
see Section lilTTI As mixing of DPPC and cholesterol molecules in 
this case is quite limited, the system probably bears traces of its ini- 
tial configuration. This applies to all state-of-the-art simulation stud- 
ies of phospholipid / cholesterol systems, and has been mentioned by 
other authors iSmondvrev and Berkowitz. 1999) . 



III. RESULTS AND DISCUSSION 
A. Equilibration 

One of the most important quantities describing lipid bilay- 
ers is the average area per molecule. The average area per 
molecule for a given configuration. A, is computed by divid- 
ing the size of the simulation box in the xy-plane, designated 
^tot, by N, the total number of molecules, i. e., DPPC's and 
cholesterols, in a monolayer. The average area per molecule 
can, among other things, be used for monitoring the equilibra- 
tion of the membrane. 

Figure [2 shows the temporal behavior of the area per 
molecule. It can be seen that after 20 ns the area per molecule 
has converged even for the highest cholesterol concentrations. 
It is, nevertheless, immediately obvious from the data that this 
type of MD simulations of bilayer systems should be at least 
of the order of tens of nanoseconds to reach equilibrium and 
surpass the longest characteristic time scales for area fluctua- 
tions. The first 20 ns of the total 100 ns were therefore consid- 
ered as equilibration, and the last 80 ns were used for analysis. 

The data clearly show the condensing effect of cholesterol: 
the area per molecule decreases with the cholesterol content. 
Further, an increasing cholesterol concentration seems to sup- 
press the fluctuations in the average area per molecule. The 
values of the average area per molecule (see also Fig.|5} are in 
excellent agreement with two recent simulation studies on the 
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FIG. 1 Temporal behavior of area per molecule. The curves cor- 
respond to, from top to bottom, cholesterol concentrations 0.0%, 
4.7 %, 12.5 %, 20.3 %, 29.7 %, and 50.0 %. 



DPPC /cholesterol system (IChiu et all l2002t iHofsaB et all 
l2003l) . As for experimental results, we are only aware of an 
accurate measurement for the average area per molecule in 
the ca se of a pure DPPC bilayer ( Naale and Tristram- 
l2000h . In this case the average area per molecule was de- 
termined to be 0.64 nm-^ at T = 323 K, in good agreement 
with (0.655 ± 0.005) nm^ obtained here. Measurements of 
the average area per m olecule in DPPC / cholesterol monolay- 
ers ( iMcConnell and Radhakrishnan, 2003) show trends simi- 
lar to ours. The exact correspondence between average areas 
per molecule measured for bilayers and monolayers, however, 
is not evident CNagle and Tristram-Nagle. 2000) . 



B. Ordering of Acyl Chains 

Average areas per mo lecule are c losely related to or- 
der parameters ( Petrach e et all Il999h . which are a mea- 
sure of the orientational order of the phospholipid tails. 
Order param eters ca n be obtained from deuterium NMR 
experi ments dSeehg a nd Seelisi 1 1 974 or computer simula- 
tions ( Tielemanetalli 1997 ). In united atom simulations such 
as ours, the orientational order can be characterized using ten- 
sors with elements Sap such that 
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Sap = 2(cos6'q cos 6*^) - -, 



(1) 



where da is the angle between t he mo lecular a axis and the 
bilayer normal (Tieleman et al.l Il997h . The molecular axes 
must be defined separately for each segment of an acyl chain: 
usually for the n^^ methylene group denoted as C„, the z axis 
points in the C„_i — C„+i direction, and C„_i, C„, and C„+i 
span the yz plane. If the motion of the segments is assumed 
to be symmetric about the bilayer normal, the experimental 
deuterium order parameter Scd can be easily acquired: 
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As the two acyl chains sn-l and s n-2 give ri se to sHghtly dif- 
ferent NMR quadrupole spHttings (iSeelig and SeeUg..l974 ). it 
is useful to compute the order parameters separately for both 
chains. 
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FIG. 2 Order parameter profiles for (a) sn~l and (b) .sn-2 tails. 
The cholesterol concentrations are: 0.0 % (o), 4.7 % (•), 12.5 % (□), 
20.3 % (■), 29.7 % (0), and 50.0 % (♦), and the index n increases 
towards the center of the bilayer. 

The order parameter profiles for the sn-l and sn-2 chains 
are depicted in Fig. |2] The ordering effect of cholesterol is 
clearly visible: the order parameters grow significantly with 
an increasing cholesterol content. For pure DPPC and low 
cholesterol concentrations, the order parameter profiles show 
a plateau for small and intermediate values of n and decay 
near the center of the bilayer When the cholesterol content 
increases, the plateau disappears, and there is a clear maxi- 
mum at intermediate n. The ordering effect of cholesterol is 
most pronounced for n ^ 6—10 and quite modest for seg- 
ments near the phospholipid headgroups and bilayer center. 
This is due to the position of the c holesterol ring system in the 
bilaye r along the bilayer normal dSmondvrev and Berkowit4 
I1999I) : the largest ordering occurs for segments at roughly the 
same depth as the ring system. For instance with 30 % choles- 
terol the order parameters for 71 ^ 6— 10 are increased roughly 
by a factor of two. 

Our results for the order pa rameters are in good agreement 
with other simulation studies ( Chiu et al., 2002: HofsaB et all 
l2003t ISmondvrev and BerkowitA .1999.) . However, as most 



force fields yield qualitatively similar results, and various 
technical details may infl uence the de tailed form of the or- 
der parameter profile ( Pat ra et alil2003l) . it is more interesting 
to make comparisons to experimental findings. 

The results for the pur e DPPC system are in good agree- 
ment with experiments ([Brown et al.L Il979t IPo uliez et aQ 
I1995L iPetrache et all 12000) As for mixtures of DPPC 
and cholesterol, Sankaram and Thompson found that when 
50 % of the DPPC molecules were substituted by choles- 
terols in a pure DPPC bilayer at T = 325 K, the order 
param eter fo r intermedi ate n was incr eased by a factor of 
2.65 dSankaram and Thompsonl ll990b). Similarly, when 
30 % of the DMPCs were replaced by cholesterols in a pure 
DMPC bilayer at T = 308 K, the order parameter increased 
by a factor of two. Vist and Davis, in turn, observed an in- 
crease by a factor of two when replacing 24 % of the DPPC 
molec ules by cholesterol at T = 323 K dVist and Davisl 
119901) . Similar agreement is found wh en our results 
are co mpared to ot her experiments ( Do uliez et al 1 119961 
Kintan ar et al 1 119861) . In all, our simulations agree well with 
experimental findings. The only detail which our, or any other 
united atom MD simulations cannot reproduce is the behav- 
ior of th e experimental deuterium order parameter for sn-2 at 
n = 2 dSankaram and Thompsonl IT990bt ISeehg and Seelid 



C. Electron Density Profiles 

Additional information about the structure of the bilayer 
along the normal or z direction can be obtained by comput- 
ing density profiles for the whole system, different molecular 
species, or certain atomic groups of interest. In simulations it 
is possible to calculate atom density, mass density, and elec- 
tron density profiles. These give information on the distribu- 
tion of atoms in the normal direction. Related information can 
be acquired from X-ray and neutron diffraction studies. Due 
to fluctuations. X-ray diffraction studies on fully hydrated bi- 
layers in a fluid phase only yield total electron density pro- 
files, whose maxi ma are associated with the elect ron dense 
phosphate groups (^N agle and Tristram-Naglell2000 ). The dis- 
tance between the maxima allows one to estimate the distance 
between the headgroups in the opposite leaflets, but does not 
yield accurate predictions for the hydrocarbon thickness or the 
true phosphate-phosphate distance (Nagle et al., 1996). Addi- 
tional information, most importantly about the average loca- 
tion of various atomic groups, can be gained from neutron 
diffraction studies either with sele ctive deuteration or in com- 
bination with X-ray diffraction llNagle and Tristram-Naglel 

iooo). 

Figure|3lshows the total electron densities calculated for the 
different cholesterol concentrations. The density profiles have 
a characteristic shape reminiscent of X-ray diffraction stud- 
ies, with maxima approximately corresponding to the location 
of the phosphate groups, and a minimum, a so-called methyl 
trough, in the bilayer center, where the terminal methyl groups 
reside. For pure DPPC and low cholesterol concentrations, the 
densities decrease monotonically from the maxima to the min- 



5 




z [nm] 

FIG. 3 Total electron density profiles as functions of distance z from 
bilayer center. The curves correspond to the various cholesterol con- 
centrations as follows: 0.0 % (dash-dotted grey), 4.7 % (solid black), 
12.5% (solid grey), 20.3% (dashed black), 29.7% (dashed grey), 
and 50.0 % (dash-dotted black). 



imum in the bilayer center. This medium density region corre- 
sponds to the methylene groups in the DPPC tails. When more 
cholesterol is present, the headgroup-headgroup distance in- 
creases, i. e., the bilayer gets thicker, and the electron density 
in the bilayer center decreases slightly. In addition, the den- 
sity in the tail region increases, and the density profile between 
the center and the headgroups is no longer monotonically de- 
creasing. The elevation is due to the fact that the cholesterol 
ring structure, which resides in the phospholipid tail region, 
has a higher electron density than do phospholipid tails. 




-3-2-10123-3-2-10123 
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FIG. 4 Electron density profiles for molecular species and atomic 
groups: (a) DPPC, (b) cholesterol, (c) water, (d) DPPC tails, (e) 
phosphorus atom (P), and (f) cholesterol ring system. The curves 
correspond to the cholesterol concentrations as indicated in Fig.|3| 



To gain more insight into the structure of the bilayer, we 
can investigate the electron densities for DPPC molecules, 
cholesterols, water molecules, phospholipid tails, phosphorus 
atoms, and cholesterol rings, portrayed in Fig.|3 All density 
profiles are consistent with a thickening of the bilayer with 
an increasing amount of cholesterol: the molecules and their 
constituent atomic groups are pushed towards the water phase. 
Still, it is clear that for all cholesterol concentrations, DPPC 
molecules largely stay within a distance of 3 nm from the cen- 
ter, while cholesterols and DPPC tails can be found within 
about 2 nm. We can conclude that cholesterol is located in 
the hydrophobic interior of the bilayer. The penetration of 
water into the bilayer becomes more difficult with increasing 
amounts of cholesterol: this reflects both the thickening of the 
bilayer, and the increasing densities in the headgroup region. 
The lipid / water interface also seems to become steeper The 
electron density of DPPC in the hydrophobic tail region de- 
creases with the cholesterol content, which is compensated by 
an increasing cholesterol electron density. By comparing the 
electron densities for cholesterol and cholesterol ring systems, 
we can conclude that only the short acyl chain of cholesterol 
can approach the bilayer center. 

Both the total electron density profile and the densities for 
molecular species and atomic groups can be compared to pre- 
vious simulations. Here we will concentrate on simulations 
on DPPC widi cholesterol at T = 323K (Hofsafi et al., 2003] 
ISmondvrev and Berkowit3ll999tlTuetal.>. .1998). In afl sim- 
ulation studies, the peaks that indicate the location of head- 
groups for pure DPPC are located approximately at the same 
distance from the bilayer center With 10- 12.5 % cholesterol 
only minor changes in the total densities can be observed. Ex- 
cept in the case of Tu et al., increasing amounts of cholesterol 
lead to a larger bilayer thickness and a slightly decreased total 
density in the bilayer center. All studies show an increased 
density in the phospholipid tail region. By investigating the 
DPPC or water densities, one may also note that all studies 
clearly indicate that the lipid /water interface becomes more 
abrupt. Our findings for the distribution of phosphorus atoms 
agree well with those of Smondyrev et al.: when the choles- 
terol content increases, the peaks are narrowed and shifted to- 
wards the water phase. In all, density profiles computed using 
slightly different force fields are, for the most part, consistent 
with each other. 

Our results are also consistent with diffraction experiments 
on DPPC and DMPC bilayers. Nagle et al. have determined 
the structure of a fully hydrated pure DPPC bil ayer in the 
liquid disordered phase using X-ray diffraction ( Nagle et all 
Il996h . The form of the density profile from our simulations 
of pure DPPC closely resembles Nagle's electron density pro- 
file for pure DPPC at T = 323 K. The head-head distance 
obtained from Nagle's experiment and that determined from 
our density profiles also are in good agreement. As for the in- 
fluence of cholesterol, Mcintosh has published X-ray diffrac- 
tion experiments on model membranes containing cholesterol 
and phospholipids with saturated tails containing 12-18 car- 
bons (Mcintosh, 1978). His DLPC/ cholesterol systems in the 
fluid phase behave in a qualitatively similar way as do our 
DPPC / cholesterol bilayers. By comparing the electron densi- 
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ties from systems with different phospholipids and cholesterol 
to the densities from pure phospholipid bilayers, Mcintosh 
also establishes the location of the cholesterol ring structure in 
the bilayer Our studies support his view. In addition, there are 
more recent neutron diffraction studies of DMPC / cholesterol 
bilayers. The studies by Douliez et al. and Leonard et al. 
clearly show that substituting 30 % of the phospholipids by 
cholesterol in a pure DMPC bilayer in the liquid disordered 
phase increases the bilayer thickness (iDouliez et ali Il996t 
lT,eonardetal. l 12001!) . Leonard et al. have also investigated 
the location of cholesterol in the bilayer, and concluded that 
cholesterol is located well within the hydrophobic core. Al- 
though DPPC has longer hydrocarbon tails than DMPC, the 
cholesterol ring structure shoul d be located in the same region 
of the bilayer (Mcintosh, 1978). Our simulations indicate that 
cholesterol is indeed situated in the non-polar region, as is the 
case in Douliez's and Mcintosh's experiments. 



D. Radial Distribution Functions 

Together, the above results ascertain that our model cor- 
rectly describes the behavior of the dimensions of the bi- 
layer and the ordering of the non-polar phospholipid tails as 
functions of the cholesterol content. Further, the structure of 
our DPPC / cholesterol bilayer in the normal direction is con- 
sistent with results from previous computations and experi- 
ments. This is very satisfactory, but in addition, we need to 
ensure that our bilayers truly are in the fluid state, i. e., that 
there is no translational long-range order This can be ascer- 
tained by examining the radial distribution functions (RDFs) 
for e. g. phosphorus and nitrogen atoms in the DPPC head- 
groups. For instance the N-N radial distribution functions 
calculated in two dimensions for various cholesterol concen- 
trations have large nearest-neighbor peaks at r 0.82 nm and 
show essentially no structure beyond r ~ 1.5 nm (data not 
shown). Additional calculations for other pairs of atoms and 
for the center of mass positions of the DPPC and cholesterol 
molecules lead to a similar conclusion, i.e., there is no lat- 
eral long-range structure. Hence, we can be confident that our 
bilayers are either in the liquid disordered or liquid ordered 
phase, as they should. With this, we consider our model to be 
valid. 



E. Estimating Average Areas per Molecule in 
Multi-Component Bilayers 

The average area per molecule, which is obtained by di- 
viding the total area of the bilayer by the total number of 
molecules, is a well-defined concept in one-component lipid 
bilayers. It includes both area actually occupied by a lipid, 
so-called close-packed area, and some free area. A similar 
quantity can be defined for multi-component bilayers. It is a 
useful quantity when simulation results are compared to ex- 
periments. Its interpretation, however, is less clear: different 
lipids and sterols could occupy significantly different amounts 
of area. Hence, it would be desirable to be able to estimate the 



average area occupied by each molecular species present in 
the bilayer. 
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FIG. 5 Average area per molecule as function of cholesterol concen- 
tration. The inset shows the average areas per DPPC (•) and choles- 
terol (o) computed as in a recent simulation study by HofsaB et al. 
The errors are smaller than the markers. 

The average area per molecule {A) as a function of choles- 
terol concentration x is portrayed in Fig.|5l As mentioned in 
Section lTlI.AI it is evident that {A ) decreases with the choles- 
terol content, a nd that the results agree well with previ ous sim- 
ulation shidies dChiu et alll2002tlHofsaB et al.Ll2003h . 

We would not, however, like Chiu et al., dare to hazard 
a guess that (A) decreases linearly with x ^nd use this as- 
sumption to compute the average areas per phospholipid and 
cholesterol. It is not obvious, in the first place, that the aver- 
age area per cholesterol or DPPC is independent of cholesterol 
content, as these authors seem to imply. 

Another way to divide the total area betwe en DPPC and 
cholesterol molecules has also been suggested (iHofsiiB et all 
'2003'). By computing the total area and volume of the simu- 
lation box as functions of the cholesterol content and making 
a number of assumptions, one can arrive at estimates for the 
average areas occupied by DPPC and cholesterol molecules. 
In this case, an important assumption is that the average vol- 
ume of a cholesterol molecule can be, for all concentrations, 
taken to be the volume occupied by a cholesterol molecule in 
a cholesterol crystal. Further, it is assumed that all space is 
occupied by DPPC, cholesterol, or water, i. e., that there is no 
free volume or area. The average areas per DPPC and choles- 
terol, aj^ppc ™d a|^Qj, obtained along these lines from our 
data, are shown in the inset of Fig.|5] These closely resemble 
the corresponding results by HofsaB et al. 

A yet further method of distributing the area among the 
molecular species in a bilayer is to apply Voronoi analysis 
in tw o dimensions (Patraetal., 2003; Shinoda and Okazak| 
Il998h . In Voronoi tessellation for a bilayer, the center of mass 
(CM) coordinates of the molecules comprising the bilayer are 
projected onto the xy plane. An arbitrary point in this plane 
is considered to belong to a particular Voronoi cell, if it is 
closer to the CM position associated with that cell than to any 
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other one. In this way one can calculate the total area associ- 
ated with the CM positions of e. g. the DPPC molecules and 
then scale this quantity by the number of DPPC molecules 
in a monolayer. The resulting average areas per DPPC and 
cholesterol, a^ppf-, and Cqij^j, as functions of the cholesterol 
content, are depicted in Fig.|6l 
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FIG. 6 Areas per DPPC (•) and cholesterol (o) computed using 
Voronoi tessellation. The errors for DPPC are smaller than the mark- 



These values for the areas per DPPC and cholesterol do dif- 
fer from those reported by HofsaB et al. This is due to the fact 
that the assumptions inherent to the respective methods lead to 
different ways of distributing the free area in the bilayer. Fur- 
ther, basic Voronoi analysis does not, in any way, allow one to 
take into account the close-packed sizes of the molecules. It 
may well be that the area, which from the point of view of the 
Voronoi analysis belongs to cholesterol, as a matter of fact, 
would be covered by projected coordinates of atoms from a 
DPPC molecule. 



F. Slicing Membranes 

We are now confronted by fundamental questions relevant 
to both one- and multi-component bilayer systems. How can 
we find estimates for the average close-packed cross-sectional 
areas for the molecular species present in a one-component or 
composite bilayer? Further, how can we estimate the average 
amount of free area in a membrane? 

Our approach to answer these questions bears a certain re- 
semblance to tomography. We map each configuration on a 
number of rectangular three-dimensional grids as follows. If 
a grid point lies within the van der Waals radius of an atom 
belonging to a DPPC molecule, this point is considered oc- 
cupied, and otherwise empty, on a grid keeping account of 
DPPC molecules. Grid points within van der Waals radiae 
of atoms belonging to cholesterol, in turn, will be occupied 
on a grid characterizing the cholesterol molecules. Finally, a 
grid for water molecules is constructed analogously. In the xy 
plane the grids have 100 x 100 elements. Because the system 



size fluctuates weakly, the size of an element will vary slightly 
from configuration to configuration. In the z direction, on the 
other hand, the size of the elements has been fixed to 0.1 nm, 
and we only consider grid points within 3 nm from the bilayer 
center. 
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FIG. 7 Cross-sections of bilayer with 20% cholesterol at 100 ns. 
DPPC grid elements have been coloured red, cholesterol is green, 
and water blue. The remaining area, i. e., the free area, is white. The 
panels correspond to slices at different distances z from the bilayer 
center: (a) bilayer center, (b) z ~ Inm, (c) z ~ 1.7 nm, (d) z ~ 
2 nm. 



The grids can be used to view given slices of the bilay- 
ers: they show cross-sections of DPPC, cholesterol, and water 
molecules, as well as patches of free area. Pictures of slices 
can be illustrative as such, and Fig. contains a selection of 
such slices for the case of 20 % cholesterol. From Figs. |3 
and 0a) we can conclude that there is quite large amounts 
of free area in the bilayer center, and that cholesterol tails 
from a given monolayer extend to the opposite monolayer. 
Panel (b) in Fig.Qportrays the region where DPPC tails and 
cholesterol ring structures should, according to Fig.|3 domi- 
nate. DPPC tails can be recognized as circular red structures, 
and the green formations are cross-sections of cholesterol ring 
structures. Part (c) is a cross-section of the bilayer at a dis- 
tance z ^ 1.7 nm from the bilayer center. Some cholesterol 
is still present in this slice, and there are also small amounts 
of water. The amount of free area is significantly smaller than 
in the bilayer center. Panel (d) finally shows a cross section at 
z 2 nm: there are DPPC headgroups, substantial amounts 
of water, and very little cholesterol. 

From the grids constructed for DPPC, cholesterol, and wa- 
ter, we can compute total area profiles for the various molec- 
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ular species, that is, average total areas occupied by the 
molecules as functions of the distance from the bilayer cen- 
ter. In addition, we can calculate free area profiles, i.e., the 
amount of free area as a function of the distance from the bi- 
layer center. In practice, this is achieved by traversing the 
grids sUce by slice and augmenting the various area profiles. 
If a grid element in a certain slice at a distance z from the 
bilayer center is occupied in, say, the DPPC grid, but not in 
the cholesterol or water grids, we increment the total area of 
DPPC in that slice, /1dppc(z)> by an area corresponding to a 
grid element. If, on the other hand, a grid point at a distance z 
from the center is occupied by neither DPPC nor cholesterol, 
nor water, the total free area Ai^^^iz) in the slice in question is 
incremented. In the end we average over the total area profiles 
constructed separately for each configuration. This procedure 
leads to a definition of free area which is similar in nature to 
the concept of empty free volume intr oduced by Marrink e t 
al. to characterize a pure DPPC bilayer ( Ma rrink et allll996l) . 
Figure UlLFl exemplifies the computation of the various area 
profiles for a bilayer with 20 % cholesterol. 



ner described above. 
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FIG. 8 Area profiles for bilayer with 20 % cholesterol scaled by 
total bilayer area: DPPC area profile {ADPPc(2))/(^tot) (solid 
black), cholesterol area profile (Achoi(2))/(^tot) (solid grey), wa- 
ter area profile {Awatcr(^:))/ (^tot) (dashed black), free area profile 
{^frcc(2:)) / (^tot) (dashed grey). The errors of the scaled areas are 
of the order of a few per cent. 



G. Close-Packed Areas for DPPC and Cholesterol 

To gain understanding of the effect of cholesterol on the 
properties of phospholipid bilayers, we first concentrate on 
the behavior of the cross-sectional area profiles for DPPC and 
cholesterol. Hence, we need to know both the total areas of 
DPPC and cholesterol and the average numbers of respec- 
tive molecules as functions of the distance from the bilayer 
center The total areas occupied by DPPC and cholesterol 
molecules, denoted by (^dppc(^)) ™d (Achoi(^)), for the 
different cholesterol concentrations are computed in the man- 
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FIG. 9 Numbers of (a) DPPC and (b) cholesterol molecules as func- 
tions of distance from bilayer center. The curves correspond to the 
cholesterol concentrations as indicated in Fig.|3| The errors are of 
the order less than 1 %. 

To find the average numbers of DPPC and cholesterol 
molecules in each slice, we locate the maximum and mini- 
mum z coordinates of each molecule with respect to the bi- 
layer center, taking into account the finite size of the con- 
stituent atoms. The molecule is considered to be present in 
all the slices between these points. By averaging over all 
molecules of a certain species and over all configurations, we 
arrive at the average numbers of DPPC molecules and choles- 
terols as functions of the distance from the bilayer center, 
denoted by (A^dppc(-2)) and (A''choi(^)), shown in Fig. |9] 
Perhaps the most notable feature in Fig. |9]is that all curves 
peak in the bilayer center. This is due to so-called inter- 
digitation: a substantial part of both DPPC and cholesterol 
molecules extend to the opposite monolayer. On both sides of 
the peak, there are broad plateaus, which reflect the amount 
of molecules of a certain species in a monolayer Eventually, 
at about 3 nm from the bilayer center for DPPC and 2 nm for 
cholesterol, the curves decay to zero. 

There seem to be two effects that together contribute to the 
thickening of the bilayer, both visible in Fig. |9] First, the 
DPPC molecules are extended. Cholesterol molecules, on the 
other hand, are not significantly elongated. These observa- 
tions are quite plausible, as the presence of cholesterol leads 
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to a smaller amount of gauche defects in the acyl chains of 
the DPPC molecules ( iHofsaB et al.Ll2003l) . Further, the tilt of 
the DPP C molecules with respect to the bilay er normal de- 
creases ( Rog and Pa senkiewicz-Gierula. '2ooiy Cholesterol, 
on the other hand, with its rigid ring structure and short tail, 
does not undergo such significant extension. 

The second effect partially responsible for the thickening 
is that with a larger cholesterol content x, a smaller amount 
of DPPC and cholesterol molecules extend to the opposite 
monolayer. Further, the ones that protrude do not penetrate 
quite as deep into the opposite leaflet as they do at low choles- 
terol concentrations. In a pure DPPC bilayer, 53 % of the 
DPPC molecules protrude to the opposite monolayer, while 
at 29.7% cholesterol the corresponding figure is 40%. The 
effect is stronger for cholesterol: At 4.7% and 29.7% con- 
centrations, respectively, 41 % and 17 % of the molecules ex- 
tend to the opposite bilayer. As the cholesterol hydroxyl is 
thought to be anchored to the DPPC headgroup via direct hy- 
drogen bonding or throu gh water bridges ( C hiu et al. . 2QD2; 
tPasenkiewicz-Gierula et a l.. 2000), this effect may be coupled 
to the elongation of the DPPC molecules. 

Equipped with the total areas occupied by the molec- 
ular species together with the average numbers of these 
molecules as functions of distance from the bilayer cen- 
ter, we can now compute the average cross-sectional 
areas for DPPC and cholesterol across a membrane, 
aDPPc(-z) = (^DPPc(^))/(^^DPPc(^)) and achoi(-z) = 
{Achoi{z))/{N^hoi{z)), shown in Fig.[Tni 

It should not come as a surprise that the cross-sectional 
close-packed area occupied by a DPPC or cholesterol 
molecule is not constant along the bilayer normal. In the 
case of DPPC, there are significant changes in the form of 
aDPPc(^) when the cholesterol content is increased. A max- 
imum located at approximately 1 nm from the bilayer center 
for a pure DPPC bilayer becomes at intermediate cholesterol 
concentrations a plateau at 0.5 — 1.5 nm from the center, and 
finally with 29.7 % cholesterol in the bilayer develops into two 
small maxima at 0.5 nm and 2 nm with a shallow minimum in 
between. 

These changes in flDPPC (z) are for the most part due to the 
behavior of the phospholipid tails: they occur in regions where 
the tail densities are high and where there are little or no head- 
groups. This can be deduced by comparing the electron densi- 
ties for DPPC molecules and DPPC tails in Figs.|4ja) and (d). 
In addition, we should note that the electron density profiles 
for the DPPC tails and the flDPPcC^) curves have many fea- 
tures common. This allows us to, at least partially, interpret 
the behavior of auppc{z) from the point of view of order- 
ing. The most substantial ordering effect with large amounts 
of cholesterol present in the bilayer occurs for carbons in the 
middle of the tail, see Fig. |2l Close to the headgroups and 
in the bilayer center the ordering effects of cholesterol are 
more modest. As increased order correlates with a decreasing 
area occupied by the tails, one expects that with an increased 
cholesterol content the cross-sectional area per DPPC approx- 
imately at a distance 1 nm from the bilayer should decrease. 
Our findings are consistent with this picture. 

This is not to say, however, that there would exist a simple 
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FIG. 10 Cross-sectional close-packed areas for (a) DPPC and (b) 
cholesterol molecules as functions of distance from bilayer center. 
The curves correspond to the cholesterol concentrations as indicated 
in Fig.0 The errors are of the order of a few percent. In the water 
phase, the relative errors for achoi are somewhat larger. 



way of mapping aoppcl^) with order parameter profiles (see 
also Section lTlI.KI i. The maximum that develops at z ~ 2nm, 
e.g., is a result of contributions from both tails and head- 
groups. From separate cross-sectional area profiles for the 
two tails and the headgroup of a DPPC molecule (data not 
shown), we found that when the cholesterol concentration in- 
creases, the cross-sectional area occupied by the tail portion 
of a DPPC molecule decreases as a consequence of ordering, 
while the area occupied by the headgroup seems to be increas- 
ing slightly (data not shown). The maximum at z ^ 2nm at 
intermediate and high cholesterol concentrations is hence re- 
lated to the interplay of the decreasing tail contribution with a 
plateau centered at z ^ 1 nm and the slightly increasing head 
contribution that peaks at z ^ 2 nm. 

In the case of cholesterol the cross-sectional close-packed 
area of a molecule is changed only weakly when the choles- 
terol concentration x is increased. The slight decrease with 
an increasing x can be explained by the tilt of the cholesterol 
molecules. At high concentrations, almost all cholesterols are 
oriented nearly parallel to the bilayer normal (data not shown). 
At low concentrations, on the other hand, the distribution of 
the angle between the bilayer normal and the ring structure be- 
comes more broad and flat, i. e., the molecules are more tilted 
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with respect to the bilayer normal. Hence the cross-sections 
appear larger at low concentrations. 

The general form of acho\{z) is compatible with our idea of 
the structure of the cholesterol molecule: narrow in the bilayer 
center where the small cholesterol tails reside and broad where 
the ring structure is located. It also reflects the thickening of 
the bilayer, as the maxima associated with the ring structures 
are pushed towards the water phase when more cholesterol 
is present. This picture, overall, supports the common belief 
that the average area per cholesterol in a phospholipid bilayer 
is largely unaltered by the amount of cholesterol in the bilayer. 

Our results for Qn hoi (z) can be compared to the o utcome of 
an old experiment iRothman and Engelmanl Il972h . where a 
model of cholesterol made of plastic was immersed in a tube 
filled with water. This experiment resulted in a steric profile 
for cholesterol, i. e., a profile of the cross-sectional area oc- 
cupied by cholesterol. This steric profile and our achoi(z), 
especially at high cholesterol concentrations, bear a surpris- 
ingly good resemblance to each other The steric profile mea- 
sured by Rothman and Engelman displays a plateau where 
the cholesterol rings are located, with cross-sectional areas of 
the order of 0.25 nm^. In the region where the cholesterol 
tail is located, they report a small maximum: here the cross- 
sectional areas are of the order of 0.15 nm^. 

It is clearly difficult to describe the close-packed area of a 
DPPC or cholesterol molecule by a single number. Of course, 
we could attempt to define the close-packed area of e. g. a 
DPPC molecule in a given DPPC / cholesterol bilayer as the 
maximum of the relevant aDPPc(2) profile, but this would 
not give accurate information about the packing of DPPC and 
cholesterol molecules in a composite bilayer. Despite this, we 
may note that the maximum values are useful at least when 
assessing the plausibility of the close-packed area profiles for 
DPPC and cholesterol molecules. 

In the case of DPPC the maxima assume values between 
0.36 nm^ and 0.42 nm^. These values can be compared to the 
average area per molecule in a pure DPPC bilayer in the gel 
state, where the contribution of the free area to the total area 
assigned to a phospholipid molecule is expected to be rather 
minor. Experiments h ave yield ed an area per mol ecule o f ap- 
proximately 0.48 nm^ ( iNagle and Tristram-Nagl j.l2000h. and 
MD s imulations suggest that {A ) — 0.46 nm^ ( Venab le et all 
l2000h . An exact comparison is not meaningful, since 
DPPC / cholesterol mixtures, especially with high cholesterol 
concentrations, have structures quite different from a pure 
DPPC bilayer in the gel state. However, the comparison 
shows that the magnitude of the close-packed areas for DPPC 
molecules is rational. 

In a similar fashion, the maxima of the acho\{z) profiles 
can be compared to values extracted from experiments on 
cholesterol crystals. The maxima found in this study de- 
crease monotonically from 0.33 nm^ to 0.29 nm^ when the 
cholesterol concentration changes from 4.7% to 29.7%. In 
a cholesterol crystal, the area per cholesterol, which in this 
case contains b oth occupied and free area, has been reported 
to be 0.38 nm^ (IChiu et al.Ll2002tlCravenIll979l:lHofsaB et all 
I2OO3I and references therein). 



H. Free Area 
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FIG. 1 1 Free areas per molecule as functions of distance from bi- 
layer center for different cholesterol concentrations. The curves cor- 
respond to the cholesterol concentrations as indicated in Fig.|3| 

We now turn our attention to the behavior of free area 
profiles for bilayers with different amounts of cholesterol. 
In Fig. ^3 we show the average amount of free area per 
molecule, i.e., afroc = {Ahcc)/N, where N is the total 
number of molecules — both phospholipids and cholesterol — 
in a monolayer. The figure clearly shows that the amount of 
free area per molecule decreases in all regions of the bilayer, 
i. e., for all values of z, with an increasing cholesterol con- 
tent. Compared to the case of pure DPPC, 4.7 % cholesterol 
in the bilayer leads to a free area per molecule reduced by ap- 
proximately 7% in all regions of the bilayer With 12.5 %, 
20.3 %, and 29.7 % cholesterol in the bilayer, the free area per 
molecule is decreased by 20 %, 35 %, and 45 %. One may 
note that the behavior of the total area of the bilayer cannot 
be explained by the reduced free area only. The occupied 
area, i. e., the area taken up by DPPC, cholesterol, or water 
molecules, also decreases with more cholesterol. For instance 
when 30 % of the DPPC molecules are substituted by choles- 
terol, the amount of occupied area decreases by approximately 
30%. 

Figure ^2 ^Iso demonstrates that an increasing cholesterol 
content in a bilayer implies that the form of the free area pro- 
file is altered. Nevertheless, the different curves correspond- 
ing to the various cholesterol concentrations have certain fea- 
tures in common: the free area profiles all have a maximum 
in the bilayer center, and there is more free area per molecule 
in the water phase than in the tail and headgroup regions. For 
pure DPPC and low cholesterol concentrations, we observe 
a minimum of free area per molecule at z 1.7 nm. For 
large cholesterol concentrations the minimum is still present, 
but due to the thickening of the bilayer it is pushed towards 
larger z: e.g. for 30% it can be found at z ~ 2nm. This 
minimum can, for all cholesterol concentrations, be associ- 
ated with a peak in the density profile of DPPC molecules 
located slightly behind the headgroups in a region where both 
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tails and headgroups are present. The density of water in this 
region is already substantial, while there is very little choles- 
terol. When the cholesterol concentration is increased, also 
another flat, plateau-like minimum starts to develop between 
the bilayer center and the minimum associated with the maxi- 
mum in the DPPC density, i. e. at z ^ 1 — 2 nm. The plateau 
is almost constant through the tail and headgroup regions. It 
has counterparts in the area profiles of DPPC and cholesterol: 
the cross-sectional DPPC area displays here a flat minimum 
and the cholesterol area a broad maximum (see Fig.llOt. We 
can thus conclude that the changes in the form of the free area 
profile are intimately related to modifications in the packing 
of the molecules in the bilayer 

It is evident that the free area profiles are related 
to the relocation and diffusion of solutes inside mem- 
branes. MD simulations suggest that solutes such 
as ub iquinone ([ Spderhall and Laaksonen, 200 ij) and ben- 
zene (Bassolino-Klimas et allll993 ) are preferentially located 
in the hydrophobic core region of a membrane. Also, it is 
known that certain non-polar probe molecules, e. g. diphenyl- 
hexatriene, prefer the bilayer center to the lipid /water inter- 
face (Lentz, 1993). These observations are in accord with our 
suggestion that the free area is largest in the bilayer center. 

There are two other simulation studies where quantities 
similar in nature to our free area profile have been calcu- 
lated for DPPC or DPPC / cholesterol bilayers. Manink et 
al. have calc ulated a so-cal l ed em pty free volume profile for 
pure DPPC dMarrink et all Il996h . This should give essen- 
tially the same information about the amount of average free 
area in a given cross-section of the bilayer as does our free 
area profile for pure DPPC. Our profile does indeed show the 
same general features as Marrink's: a maximum in the bi- 
layer center and minima near the headgroup region. Tu et 
al. have also looked at the influence of 12.5 % choleste rol on a 
so-called empty free volume fraction d Tuet al.Lfl99l . which 
is equivalent to our total free area scaled by the total area of the 
bilayer. If we compare such scaled free areas (data not shown) 
to Tu's data, we see that the scaled profiles have many features 
in common. One difference is that the bilayer thickening is not 
visible in Tu's results, while it can be clearly distinguished 
from ours. The thickening has also been verified experimen- 
tally. Further, there are some differences in the detailed form 
of the profiles in the bilayer interior, i. e., the location of the 
minima are slightly different. As Tu et al. point out, the dif- 
ferences are probably due to different computational models. 



I. Lateral Diffusion and Free Area 

We have seen that an increasing cholesterol concentration 
reduces the amount of free area per molecule in the bilayer 
and simultaneously alters the packing of the molecules. On 
the other hand, it is well known from experiments that lateral 
diffusion of both DPPC and cholesterol molecules is affected 
by changes in the cholesterol content ( Al meida et al l ll992t 
f^'ilippov et al., 2003a; Konig et al., 1992). It is reasonable to 
expect that these properties of the bilayer and the observed 
modifications in them with the cholesterol concentration are 



related. Free volume theory is a simple but appealing model 
for explaining such dependencies. 

Free volume theory was originally developed for de- 
scribi ng the transport pr op erties o f glass-forming flu- 
ids dCohen and Turnbulil Il959l [Macedo and LitovitS 
119651: iTurnbuU and Coheni 1196 iL 119701) . It was 
subsequently adapted to modeling two-dimensional 
diffusion (Almeida e t al,, 1 992; Ga ll g etal 197^ 
iMacCarthv and Kozakl Il98l IVaz et all 119851) and is 
usually in this context dubbed free area theory. Free area 
theory, a two-dimensional mean-field model for diffusion, can 
be used to at least qualitatively describe lateral self-diffusion 
in lipid bilayers dAlmeida et all Il992h . According to free 
area theory, lateral diffusion of a lipid or sterol in a bilayer 
is restricted by the occurrence of a free area greater than 
some critical area adjacent to the diffusing molecule. A 
diffusing molecule spends a compa ratively long time — of 
the order of tens of nanoseconds dTieleman et all Il997t 
Vattulainen and Mouritsen, 2003) — in a cage formed by its 
neighbors, and then, given a large enough activation energy 
and an adjacent free area, jumps. 

More specifically, free area theory predicts that the lateral 
diffusion coefficient of a lipid or sterol diffusing in a bilayer 
depends on the free area and the packing properties as fol- 
lows (Almeida et al., 1992): 

~ exp(-ao/af). (3) 

Here ao is an estimate for the average cross-sectional area for 
a DPPC or cholesterol molecule and af is a measure for the 
average amount of free area per molecule in the bilayer 

To examine the validity of Eq. (|3} we compute the lateral 
diffusion coefficients for DPPC and cholesterol molecules at 
different cholesterol concentrations. The lateral tracer diffu- 
sion coefficients can be computed using the following Einstein 
relation 

- l™ JT^r E - ^^'(0)]')- (4) 

t^oo ^A/spccios ^ 

^ 2 — 1 

Here n (t) is the CM position of molecule i at time t and the 
sum is over all molecules of a given species. The lateral dif- 
fusion coefficients have been calculated by following the po- 
sition of each molecule in the upper (lower) monolayer with 
respect to the center of mass position of the corresponding 
upper (lower) monolayer Thus, should there be any drift, the 
motion of the CM of each monolayer has been taken into ac- 
count. 

Results for lateral diffusion coefficients are shown in 
Fig. The lateral diffusion coefficients for both DPPC 
and cholesterol decrease monotonically with an increasing 
cholesterol content. T his reduction i s qualitatively consis- 
tent with experiments ("Almei da et all 1 19921 'Filippov et all 
I2OO3 a; Konig et al., 1992). Quantitative comparisons should 
preferably be made to experimental techniques that probe lat- 
eral diffusion of individual molecules at time scales compara- 
ble to those reached in MD simulations. Fluorescence correla- 
tion spectroscopy (ECS) measurements should hence give us a 
good reference. In ECS measurements for DLPC/ cholesterol 
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FIG. 12 Lateral diffusion coefficients of DPPC (•) and cholesterol 
(o) molecules as functions of cholesterol concentration. 



systems Korlach et al. found that when the cholesterol con- 
centration was increased from 0% to 60%., I?t for DLPC 
was reduced by a factor of ten (HCorlachet al., 1999). Even 
though the acyl chains of DLPC molecules are shorter than 
those of DPPC molecules, our findings are in reasonable ac- 
cord with Korlach's experiments. 

Let us now consider the implications of our results to free 
area theory for lateral diffusion. In free area theory, the criti- 
cal area ap is essentially a number describing the close-packed 
cross-sectional molecular area of the diffusant. In the same 
spirit, the average free area per molecule Of should be charac- 
terized by a single number. We have, however, seen that the 
free areas per molecule and the areas per DPPC and choles- 
terol molecules are functions of the distance from the bilayer 
center Hence, it seems that a two-dimensional mean-field 
model might be a too simplistic means of describing lateral 
diffusion. 

In our opinion, one should at least not expect free area 
theory to yield quantitative results. It might, however, give 
qualitative predictions about trends in cases where e. g. the 
cholesterol content in a bilayer is increased. With this in 
mind, let us assume that the cholesterol concentration in a 
DPPC /cholesterol bilayer rises from 4.7 % to 29.7 %. If we 
now use the largest possible values for the fractions ag/af, 
free area theory will give us upper bounds for the reduction 
of the lateral diffusion coefficients. The lateral diffusion co- 
efficient for DPPC should according to free area theory be 
now reduced by a factor three at most, and I?t for cholesterol 
should decrease by a factor of two. As a matter of fact, the 
lateral diffusion coefficients for both DPPC and cholesterol 
computed from the simulation data are reduced much more 
strongly, see Fig. [2] We can conclude that Eq. (|3} tends to 
underestimate the changes in the values of the lateral diffusion 
coefficients. 

Even though the discrepancies in the predictions of Eq. (|3} 
and the computed lateral diffusion coefficients do exist, we 
cannot immediately declare free area theory incomplete. 
There is a detail that has been overlooked in our discussion so 



far, and the significance of this detail will now be considered. 
To jump to an adjacent empty site, a diffusing molecule needs 
energy to overcome an activation barrier. In free area theory 
this is accounted for by letting the lateral diffusion coefficient 
be proportional to a Boltzmann factor exp(— £'a/ k-oT), where 
i?a is the activation barrier. As a growing cholesterol concen- 
tration increases the ordering of the DPPC tails and therefore 
reduces the area per molecule, it seems reasonable to expect 
that i?a should increase with the cho lesterol content. Experi- 
mental results (! Almeida et aT*. '19921 do support this idea but 
are partly contradicting. This is, however, prob ably due to 
the fitting procedure used ( Almeid a et al l ll992h . In a more 
recent study, Filippov et al. used NMR to study the lateral dif- 
fusion in palmitoyloleoylphosphocholine/ cholesterol and di- 
oleoylphosphocholine/ cholesterol bU ayers over a cholesterol 
concentration range of '^O - 45 mol % dFihppov et al., 2003} 3j). 
At small X, they found the apparent (Arrhenius) diffusion bar- 
rier to be approximately constant, while for large x the dif- 
fusion barrier increased markedly. Hence, the neglect of the 
energy term might in our case lead to slight underestimates for 
the reduction of the lateral diffusion coefficients. 

Summarizing, we have found that free area theory correctly 
predicts the reduction of the lateral diffusion coefficients with 
an increasing cholesterol concentration. At the same time it 
seems unnecessary to aim for a quantitative description with 
such a simple framework. Instead of being based on mean- 
field arguments, a full theoretical description of lateral diffu- 
sion should account for local free volume fluctuations in the 
vicinity of diffusing molecules. Atomic-scale MD studies in 
this direction should be feasible in the near future. 



J. Area Compressibility 

Lateral diffusion is clearly influenced by the average 
amount of free area in the bilayer. However, not only the aver- 
age free area, but also fluctuations in the amount of free area 
should play a role here. Recall that free area theory states 
that a diffusion jump is not possible unless there is a larg e 
enough free area next to the diffusant ("Alme ida et all [T992i) . 
Large enough free areas are a result of fluctuations, and hence 
we would expect diffusion to depend on the magnitude of the 
fluctuations: decreasing fluctuations and slowed lateral diffu- 
sion should be coupled. For similar reasons, it is likely that 
permeation of molecules across membranes can at least par- 
tially be explained by area fluctuations in membranes. 

We may quantify area fluctuations in different regions of 
the membrane as follows. The starting point is the average 
occupied area {Aocc{z)), i. e., the area which is not free but 
occupied by DPPC, cholesterol, or water molecules. The oc- 
cupied area obviously varies with the distance from the bilayer 
center z. Based on the occupied area we define an area com- 
pressibiUty as follows: 



{Aocc{z)) 



(5) 



Here fee is the Boltzmann constant and {6A^^^) — {A^^^ — 
(Aocc)^- The area compressibility is a measure of the fluctu- 
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ations in the occupied area: a high compressibility indicates 
small fluctuations and a low compressibility, correspondingly, 
large fluctuations. Hence, the area compressibility should be 
related to the permeation of small solutes, as well as for the 
lateral diffusion of lipids and sterols. 
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FIG. 13 Area compressibilities as functions of distance from bilayer 
center. The curves correspond to the cholesterol concentrations as 
indicated in Fig.|3| The errors are between 10 and 20 %. 

Figure^lshows the area compressibility profiles computed 
for systems with different amounts of cholesterol. Before fo- 
cusing on the behavior of ka{z), let us stress that these quan- 
tities are sensitive to force fields and various computational 
details and can only be used for discussing qualitative trends 
with an increasing cholesterol content. 

Regardless of the cholesterol content, all compressibility 
profiles have a minimum in the bilayer center. Moreover, the 
values of compressibility are identical in the center. The situ- 
ation in the water phase makes sense: we expect that the com- 
pressibilities, irrespective of the cholesterol content, should 
be approximately similar The interesting regions are the tail 
and head ones. The compressibility profiles show two maxima 
between the bilayer center and the water phase, the first at ap- 
proximately 1 nm from the bilayer center and the second at 
1.7 — 2.0 nm, depending on the cholesterol concentration. Be- 
tween these we observe a local minimum. For pure DPPC and 
at low cholesterol concentrations there is a very flat, plateau- 
like maximum centered at 1 nm. With more cholesterol, the 
maximum grows considerably. Returning to Fig.|3f), we note 
that the position of the growing maximum coincides with the 
location of the cholesterol ring structure. Therefore we can 
conclude that the cholesterol steroid rings strongly reduce the 
area fluctuations in the bilayer. From the point of view of 
free area theory, the region with the ordered DPPC tails and 
cholesterol rings seems to be the rate-limiting region for lat- 
eral diffusion of lipids and sterols. 

The local minimum between the two maxima moves from 
a distance 1 .5 nm from the bilayer center to 1 .9 nm from the 
center. This means that the minimum is located in a part of the 
bilayer with the uppermost tail methylene groups and parts of 
headgroups. The densities of cholesterol backbone are quite 



small here, while the electron densities of the cholesterol hy- 
droxyl groups peak (data not shown). Very few cholesterol 
rings, and hence tails with less order than at 1 nm, and pos- 
sibly also the interface between hydrophobic and hydrophilic 
parts, lead to slightly larger area fluctuations here. This could 
have consequences for the permeation of small molecules. 
Knowing that there are larger area fluctuations in this region 
than elsewhere, does not, however, tell us how permeation is 
affected. Jedlovszky et al., e. g., found in a recent simulation 
study of DMPC / cholesterol that the region with the choles- 
terol hydroxyl groups is indeed important from th e point of 
view of permeation ( Jedl ovszkv and Mezeil 120031) . The ef- 
fect on the actual rate of the permeation process, nevertheless, 
must depend on the properties of the permeant molecule. 



K. Close-Packed Areas from Ordering of Acyl Chains? 

Let us return to the average area per DPPC and investigate 
whether anything can be said about the close-packed area of a 
DPPC molecule based on the chain order parameters. Tradi- 
tionally, the use of deuterium NMR experiments to determine 
the average area per DPPC has resulted in a wide variety of 
values (Nagle and Tris^tram-Nagle, 2000). This is not so much 
due to the underlying results for the order parameters as due 
to the interpretation of the results. Petratche et al. have rather 
recently suggested a way of relating the deuterium order pa- 
rameter to_flie_average chain travel distance along the bilayer 
normal (iPetrache et al.Lll999h : 
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(6) 



Here (I?„) is the average chain travel distance along the bi- 
layer normal for segment n, Dm the maximum travel per 
methylene for all-trans chains oriented perpendicularly to the 
bilayer, and the deuterium order parameter for segment 
n. By assuming that (A„) « Vch2 /i^n), where VcHa is 
the volume per methylene group (this is only true if (I?^) ~ 
{Dn)"^, see (Petrache et aL.1999.) ) and recalling Eq. (|2|i, we 
may write: 
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where Aq is the area occupied by a fully ordered phospholipid 
molecule. By examination of Figs.l2land ll0r a). we can con- 
clude that it is unrealistic to expect that Eq. Q should allow 
one to extract the detailed form of aDPPc(z) from S^z- Nev- 
ertheless, Eq. might be useful in predicting the average 
areas per DPPC molecule in the tail region, e. g. at a distance 
1 nm from the bilayer center, where the headgroup density is 
negligible for all cholesterol concentrations. 

To find the values of the order parameters at 1 nm from 
the bilayer center, we use electron density profiles calculated 
separately for each methylene group in the hydrocarbon tails 
(data not shown) to determine which segment is located at a 
distance 1 nm from the center for each cholesterol concentra- 
tion separately. The order parameters at 1 nm are then cal- 
culated as averages over the segments whose electron density 
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profiles peak at the close vicinity of 1 nm and over the sn-l 
and sn-2 tails. The close-packed areas for DPPC at 1 nm from 
the center, in turn, can be easily obtained from the aDPPc(^) 
profiles. The resulting values and a fit to Eq. Q are shown 
in Fig. ^5 The fit is astonishingly good, given that Eq. 
has been developed for a pure phospholipid bilayer and is 
based on a rather simple model. The best fit is obtained with 
Aq w 0.28 nm^. Aq should in this case be interpreted as the 
area occupied by the fully ordered sn-l and sn-2 tails. Hence, 
the agreement with Fig. II Uf a) is surprisingly good. Yet one 
should not pay too much attention to the exact numerical value 
here, as it probably depends on the details of the force field. 

0.6 I — I — 1 — I — I — I — 1 — 1 — I — 1 — 1 — I — 1 — I — I — 1 — 1 — I 
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FIG. 14 Order parameters vs. close-packed areas for DPPC at 1 nm 
from bilayer center. The markers represent values computed from the 
simulations and the solid line is a fit to these data based on Eq. 0. 

H ofsaB et al. have use d Eq. in a slightly different set- 
ting ( HofsaB et al.ll2003h . As order parameters they have used 
averages over the order parameter profiles from segment 3 to 
segment 8, and the average areas per DPPC they used contain 
some free area. They, too, find that Eq. gives a very good 
fit to their data. However, we expect that order parameters are 
related to close-packed cross-sectional areas for DPPC chains 
rather than to average areas per DPPC containing an arbitrary 
amount of free area. 



the bilayer. Using this method we have computed the aver- 
age cross-sectional areas for DPPC and cholesterol, as well as 
the total free area, as functions of the distance from the bi- 
layer center The method should be generally applicable for 
all kinds of pure and composite bilayers. Moreover, it could 
be used for investigating bilayers with integral proteins and in 
such a way finding out how the bilayer structure is changed in 
the vicinity of embedded proteins. 

Inspection of the cross-sectional close-packed area profiles 
for DPPC and cholesterol, i. e., the close-packed areas as func- 
tions of the distance from the bilayer center, has shown that 
cholesterol alters the packing of molecules and reduces the 
amount of occupied space. These phenomena have been quite 
generally explained in terms of the form of the cholesterol 
molecule and the ordering effect of cholesterol on parts of the 
phospholipid tails. 

Cholesterol has also been found to significantly reduce the 
average amount of free space in all regions of the bilayer. We 
have further discovered that the form of the free area profiles, 
i. e., the average amount of free area as a function of the dis- 
tance from the bilayer center, is altered. These changes seem 
to reflect the ones observed in the close-packed area profiles 
for DPPC and cholesterol. We therefore conclude that the 
packing and free area properties are strongy coupled. 

Also lateral diffusion of DPPC and cholesterol molecules 
has been found to be strongly reduced with an increasing 
cholesterol content. Further, the changes in the packing prop- 
erties and the average amount of free area seem to be reflected 
in the behavior of the lateral diffusion coefficients for DPPC 
and cholesterol molecules. We have, however, learned that 
even though so-called free area theories correctly predict the 
suppressed lateral diffusion with reduced free area, the depen- 
dence cannot be quantitatively described by mean-field mod- 
els such as free area theory. Not only are the average free 
areas or volumes relevant for diffusion, but also the size distri- 
bution, shape, and local fluctuations of the free volumes in the 
bilayer are important. It would hence be interesting to see how 
cholesterol influences the size distribution of free volumes in 
the bilayer. 



IV. SUMMARY AND CONCLUSIONS 

We have performed 100 ns molecular dynamics simula- 
tions at T = 323 K on a pure DPPC bilayer and compos- 
ite DPPC /cholesterol bilayers with 4.7%, 12.5%, 20.3%, 
29.7 %, and 50.0 % cholesterol. The main focus has been 
on the packing of molecules, free area in different parts of 
the bilayer, and lateral diffusion of DPPC and cholesterol 
molecules. Especially the interplay between these properties 
has been considered. 

To investigate the packing and free area properties, we have 
introduced a novel method for estimating the average space 
occupied by DPPC, cholesterol, and water molecules, along 
with the average amount of free space, in different regions of 
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